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ABSTRACT 



Aims. We analyzed the resolved stellar population of the C component of the extremely metal-poor dwarf galaxy I Zwl8 in order to 
evaluate its distance and star formation history as accurately as possible. In particular, we aimed at answering the question of whether 
this stellar population is young. 

Methods. We developed a probabilistic approach to analyzing high-quality photometric data obtained with the Advanced Camera for 
Surveys of the Hubble Space Telescope. This approach gives a detailed account of the various stochastic aspects of star formation. 
We carried out two successive models of the stellar population of interest, paying attention to how our assumptions could affect the 
results. 

Results. We found a distance to the C component of I Zw 18 as high as 27 Mpc, a significantly higher value than those cited 
in previous works. The star formation history we inferred from the observational data shows various interesting features: a strong 
starburst that lasted for about 15 Myr, a more moderate one that occurred » 100 Myr ago, a continuous process of star formation 
between both starbursts, and a possible episode of low level star formation at ages over 100 Myr. The stellar population studied is 
likely « 125 Myr old, although ages of a few Gyr cannot be ruled out. Furthermore, nearly all the stars were formed in the last few 
hundreds of Myr. 

Key words. Galaxies: individual (I Zw 18), Galaxies: photometry, Galaxies: stellar content, Galaxies: formation 



1. Introduction 

The study of dwarf galaxies is an important topic for the un- 
derstanding of how galaxies form and evolve. Their presence in 
galaxy clusters and their physical properties put constraints on 
cosmological mode ls, especially on th e dark matter content of 
the Universe (e.g. iRobertson et al1l2005l) . Furthermore, because 
their shallow gravitational potential well, dwarf irregular (dlrr) 
galaxies are ideal benchmarks for studing the effects of various 
processes, either internal (e.g., stellar winds and supernovae) or 
external (like gas infall and tidal forces), o n the triggering and 
regulation of star formation (SF) in galaxies dHunter et al. 1119981 
iHensler & Ries chick 2002). Finally, metal-deficient dlrr galax- 
ies are expected to match the chemical composition of pristine 
galaxies, and efforts have been made to derive the primordial 
helium abundance of t he Universe from the spectroscopic analy- 
sis of their ionized g as (Peim bert & Torres-Peimbertl 1974 [19761 
llzotov et al.iri997bl 120071 iPeimbert et al.ll2007l) . 

Among the dlrr galaxies known, I Zwl8 is arguably the 
most fascinating one. This object is the second lowest metallic - 
ity galaxy known (12 + log O / H = 7.2), corresponding to » 1/30 
of the solar value (e.g . Page! et al. 1992 , ISkillman & Kennicuttl 



1993| llzotov & Thuanll 19981 llzotov et al.lTl 999). bein g the low- 
est o ne SBS 0335-052 (West) (llzotov et all 12004). Its blue 
color (Ivan Zee et aT 1998) a nd its intense nebular em ission (e.g. 
IVflchez & Iglesias-Paramoll998llCannon et al.ll2002) bears wit- 



ness to an intense ongoing SF episode. The galaxy presents a 
complex morphology; it is dominated by a two-lobed body (the 
"main body") with one companio n, the "C component" (here- 
after I Zwl8C lDufour et al.lll996ah . 

In spite of many efforts, the distance and, more important, 
the age of the galaxy have not been well determined yet. By 
analyzing HST/WFPC2 color-magnitude diagrams (CMDs) and 
HST/FOS spectra, and assuming a distance of 10 Mpc (apply- 
ing straightly the H ubble law to the redshift of the galaxy), 
iDufour et al.l (Il996bl) detected a stellar population of up to 50 
Myr in the main body of the galaxy and up to 300 Myr in the 
C component. They found no support for or against an old un- 
derlying population. Assuming the same distance and examin- 
ing a larger set of HST/WFPC2 photometric data, lAloisi etail 
(1 19991) decomposed the star formation history (SFH) in a con- 
tinuous SF lasting over ~ 1 Gyr and an ongoing starburst that 
started * 15-20 Myr ago. By analyzing HST/WFPC2 CMDs 
jointly, ground-bas ed spectra and the morphology of some neb- 
ular features, llzotov et all (1 19991) find that a distance of 20 Mpc 
is necessary to explain the ionization state of the gas. They 
also state that the age of the galaxy, as derived from all three 
kinds of data, d oes not exceed 100 M yr, in agreement with the 
conclusion from lIzotov&TTiuanl d 19991) th at very metal-deficient 
galaxies must all be young. lOstlinl (120001) carried out an analy- 
sis of HST/NICMOS photometric data and detected asymptotic 
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giant branch stars (AGB) that, at the adopted distance of 12.6 
Mpc (derived from the redshift corrected for the Virgocentric 
flow), are at least 1 Gyr old. From the surface dist r ibutio n of 
the fluxes and colors of the galaxy, iKunth & Ostlinl d2000l) ar- 
gue that I Zwl8 is an old galaxy whose age is likely ~ 5 Gyr. 
A signific ant age of po s sibly several Gyr has also been advo- 
cated by iLegrand et al.l d2000t) to account for the very homo- 
geneous distribution of heavy elements throughout the ionized, 
optical-emitting gas of the galaxy. By studying the deep spec- 
tra of I Zwl8C. llzotov et all I d2001l) find that a distance of ~ 15 
Mpc is needed to reconcile ionized gas in this object with the 
apparent magnitudes of its brightest blue stars. This distance can 
actually be interpreted as a lower limit, since at larger distances 
these stars would be still considered as ionizing. Their spectra 
are well-fitted by models with a continuous star formation ex- 
tending over the past 10-100 Myr, with no evidence o f any older 
underlying stellar population. Pap aderos et al.l ((2002) compared 
surface brightness maps obtained through various optical and 
near-infrared filters with the HST, which did not reveal any stel- 
lar population older than 0.5 Gyr. Such a population may exist, 
but most of the stellar content must be younger than this age. 
llzotov &Thuanld2004l) carried out deep HST/ ACS observations 
of I Zwl8, improving by « 2 mag the depth of stellar photo- 
metric data with respect to previous efforts. They concludes that 
the distance to I Zwl8 should lie somewhere between 12.6 and 
15 Mpc with the latter value being more likely, and claim that 
the absence of any detected red giant is proof that the galaxy is 
young (< 500 Myr). F inally, by analyzing the same images as 



Izotov & Thuan (2004), along with other deep HST/ ACS data, 



Aloisi et al.l (120071) detected a red giant branch (RGB), as well 
as a confirmed, Cepheid star. They derived a distance of » 18 
Mpc and argue that the presence of RB G stars shows th at the 
galaxy is old (at least 1 Gyr, but see also iTosi et a 020071) . 

It is striking to see that similar approaches to the evalua- 
tion of the a ge of I Zwl8 lead to contradict ory results (e.g., 
Dufour et aljll996bl vs. lAloisi et al.l 1 19991 or IKunth & O stlin 



2000 vs. Papadero s et al. 1 120021) . It turns out that the age esti- 



mates from the authors mentioned above go from a ten Myr to 
several Gyr! This might be explained at least partially by some 
uncertainties in the studies summarized above, such as the lim- 
ited depth of the available data or the sensitivity of CMD anal- 
ysis to the distance and extinction adopted. For this reason, we 
decided to carry out a ne w analysis of the deep HST/ ACS ob- 
servations performed by llzotov & Thuanl d2004l) . focusing on 
I Zwl8C. The latter has been little studied because of its faint- 
ness (« 24 mag/n" in the V band). However, there is a key dif- 
ference between I Zw 18C and the main body that makes the C 
component more attractive for CMD studies: it is older. 

1 . There is a significant absence of nebular emission. Because 
it is older, the observed surface brightness in the C compo- 
nent is lower than in the main body and is predominantly 
produced by the stellar component. The surface brightness 
in the main body is prod uced by both the stellar and nebu - 
lar components (see, e.g. IVflchez & Iglesia s^aramof l998). 
Also, compact clumps of ionized gas in the main compo- 
nent would mimic stellar sources, and the patchy ionized gas 
emission makes background subtraction more difficult. 

2. Isochrones (and their associated stellar luminosity functions) 
have a simpler structure at older (post Wolf-Rayet) ages, with 
well-defined structures excepting the blue loops that make a 
small contribution in the stellar luminosity function. In the 
case of the main component, there are identifie d Wolf-Rayet 
stars dLegrand et al.ll 19971 llzotov et all 1 1 997al) such a pres- 



ence is difficult to explain in low- metallicity environm ent 
with standard evolutionary models iBrown et al.l d2002l) . so 
more complex evolutionary models (like the inclusion of 
rotation) would be required. However, current models in- 
cluding rotation produce evol utionary tracks with com plex 
structures (see e.g. Fig. 3 in iMevnet & Maederl 120051) . so 
isochrones becomes strongly dependent on the assumed ini- 
tial rotational velocity, the possible distribution of rotational 
velocities, and the (still unexplored) validity of homology re- 
lations between different stellar tracks for rotating stars. 

There is an additional reason for choosing the C component 
instea d of the main body. In the distance obtained bv llzotov et al.l 
(2001) from the integrated spectra of the galaxy the authors take 
the IM F sampling effects in the emission line spectrum into ac- 
count dCervino et al.ll2000l in its first version). However, they 
did not c onsider IMF sampling effects in the analysis of th e con- 
~~ aD[2002, ICerv ino & Luridianal l2006) from 



tinuum dCervino et al 
which the distance of 15 Mpc to I Zwl8C (and by extension, 
the main body) was obtained. It was also our intention to obtain 
an alternative value for the distance with a method where such 
sampling effects would be minimized (the CMD analysis). 

For these reasons, we have decided to use study I Zwl8C 
instead of the main body. We had built a new, improved CMD 
from the images of Izotov & Thuan ( 2004) and analyzed it with 
a probabilistic approach to evaluate the distance, the SFH, and 
the age of this component. 

This paper is divided as follows. In Sect. [2] we define the 
probabilistic tools used in our approach. In Sect. [3] we describe 
the observational data and the theoretical isochrones on which 
our work was based. In Sect. [4] we present a first model of the 
stellar population studied. Section|5]is dedicated to a revision of 
the assumptions of this model, and in Sect. [6] we present a new 
model with improved reliability. We infer the minimum age of 
I Zwl8C and discuss the estimate of its distance. Finally, we give 
our conclusions in Sect. [7] 



2. A probabilistic approach 

In considering a stellar population in which N individual stars 
are detected in a set of n photometric bands, note m U j the mag- 
nitude of the jth star in the z'th photometric band. We first ex- 
amine the case of a simple stellar population, i.e. a population 
where all the stars were formed at the same time and in the same 
conditions. The initial mass Mj of the jth star is considered a 
random variable and the initial mass function (IMF) <f>(M) is 
regarded as the probability density function (PDF) from which 
Mj is drawn. The expected magnitudes my of the star depend 
on Mj and on properties of the stellar population: the age r, the 
metallicity Z, the distance d and the interstellar extinction coef- 
ficients A,. Furthermore, those magnitudes are affected by pho- 
tometric errors, which we suppose follow Gaussian distributions 
of widths try. Finally, the photometric data do not usually cover 
the whole IMF of the stellar population, since the faintest stars 
are not detected. Consequently, if we call p c (m) the photometric 
completeness function of the data (the probability a star of mag- 
nitudes nti is detected), then the observed magnitudes my of the 
star j are associated to the following likelihood function (LF): 



Lj(r,Z,d,A) = 
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with 



<p'(M\T,Z,d,A) = 



and 



<p(M)p c (m((M)) 



J^f <P(M')pc(m(M'))dM' 



~ -At (m.; : - mi ,) 2 

X*(M;T,Z,d,A) = Yj 



(2) 



(3) 



In Eqs. Q]and[2] the upper limit /Vt up (T) of the integrals is the 
maximum initial mass with which a r-old star still visible today 
could have been born. 

Assuming that the initial masses of the stars are mutually 
independent, the LF L pop associated to the distribution of magni- 
tudes of the ;V stars is given by the product of the individual LFs 



Lp 0p (r, Z, d, A) = Y\ Lj(r, Z, d, A). 

7=1 



(4) 



In practice, it is preferable to work not with the LFs themselves, 
but rather with their logarithms. Hence, we introduce the follow- 
ing quantities: 



"Hpop = log Lp p = ^ r H j . 



(5) 
(6) 



In this section, we have limited the definition of Lj to the 
case of a simple stellar population. In what follows, we propose 
another definition that can be applied to composite populations. 



2.1. Introducing the stellar formation history into the entropy 

Many stellar populations show prolonged periods of star forma- 
tion, in which case Eq.[T]cannotbe used. It is necessary to intro- 
duce a function N(t) that describes the SFH of the population 
studied. Like their initial masses, the ages of the different stars 
are supposed to be independent random variables. As a result, we 
define N(t) dr as the expectancy of the number of stars observ- 
able today and whose ages are comprised between t and t + dr. 
Consequently, Lj(T,Z,d,A) must be replaced by the following 
LF: 



L'j(N;Z,d,A) = 



Ij(N;Z,d,A) 

j N(T)dT 



with 



Ij(N;Z,d,A) = j N(t) Lj(t; Z, d, A) dr. 
The entropy fipop hence become: 



pop 



n „ 

Y_}ogIj(N;Z,d,A)- N\og I N(j) At. 
;'=i ^ 



(7) 



(8) 



(9) 



For given values of Z, d, and A, maximizing < K pop 



leads to 

the best-fit shape of N(t), not to absolute values of this func- 
tion. The latter must then be normalized with the number of stars 
used: 



N(t) dr = N. 



(10) 



The mass rate A1(t) of star formation can be derived from N{t) 
by considering which part of the IMF is visible in the data at the 
age r: 



M( T ) = N(t)- 



rAl„p(0) 



M<f>(M)dM 



f™f } 4>(M)Pc(m(M;r))dM 



(ID 



In some cases, a parametric model can be provided for 
the SFH. In this case, 7V(t) can be fitted using a Levenberg- 
Marquardt algorithm. In the case where no parametric model 
is assumed for the SFH, the latter must be inferred in a self- 
consistent way from the observed distribution of magnitudes 
m, ;. This can be ach ieved with the method developed by 
iHernandez etalld 19991) . 

2.2. Comparison with other approaches 

The variable "T/pop is an entropic tool whose maximization can 
be used to analyze resolved stellar populations. Contrary to sim- 
pler estimators (e.g., the chi-square one), it is sensitive not only 
to the shape of the tested isochrones in color-magnitude dia- 
grams, but also to the density distribution of the stars along 
those isochrones, through the relations between the age, mass 
and magnitudes. This is particularly true for rapid stellar phases, 
where the magnitudes are very sensitive to the ages and initial 
masses of the stars. Other probabilistic methods have been de- 
veloped to characterize resolved stellar populations, especially 
the comparison of star counts between observed and synthetic 
Monte-Carl o CMDs in differe nt bins of the color-magnitude 
space (e.g. iGallart et al.l 1 1 9991) . However, the formalism pro- 
posed in this work is based on rigorous concepts, it is more ac- 
curate than the synthetic CMD approach. In particular, it avoids 
the biases and errors caused by the binning of data and the ran- 
dom content of simulated CMDs. Moreover, it bypasses expen- 
sive Monte-Carlo computations for the search of optimized pa- 
rameters of the stellar population studied. In summary, < 7/ pop al- 
lows a study of the stellar luminosity function as a whole, with 
its correlations among different luminosity bins, instead of only 
portions of the luminosity function without correlations among 
bins. 

The approach we have presented is of course not free of 
drawbacks. In particular, there is no simple way to assess the 
goodness-of-fit, and Monte-Carlo simulations may be required 
to carry out this task (see Sect. 14.31 ). Furthermore, like other 
probabilistic methods, the results deeply depend on the input 
physics, in particular the stellar evolution models and the conver- 
sion of stellar parameters to magnitudes. Ideally, different evo- 
lutionary tracks and model atmospheres, or empirical magnitude 
calibrations, should be tested when analyzing an object, at the 
cost of simplicity and time. 



3. Applications to I Zw18C 

3.1. The photometric data 

We carried out photometric measurements of the HST/ACS im- 
ag es obtained in the bands F555W (~ V; Fig.ffl) and F814W (~I) 
by llzotov & Thuanl d2004h . They performed photometric anal- 
ysis of those images, obtaining data significantly deeper than 
those of previous works. However, their procedure, made use 
of a model point spread function (PSF) that was not optimal. 
For this reason, we performed our own measu rements by us - 
ing the well-suited DOLPHOT/ACS package dDolphinl l20Q0). 
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Fig. 1. HST/ACS image of I Zwl8C through t he F555W filter 
The tw o clusters "C" and "NW" mentioned bv llzotov & ThuarJ 
(2004) are reported. The image display scale is logarithmic. 
North is up and East to the left. The "white pixels" are contami- 
nated pixels that were rejected when performing the photometric 
measurements. 



Fig. 2. Color-magnitude diagram of I Zwl8C. The typical error 
bars as a function of F814W are reported for both the blue and 
red branches. The two clusters are annotated with open circles. 
The stars marked with crosses are those we rejected from the 
useful sample due to their extreme red color. We also report the 
90% and 50% completeness limits as dashed lines. 



which correctly accounts for the pixel undersampling of the in- 
struments and uses accurate model PSFs. 

The CMD that we obtained is shown in Fig. [2] The F555W 
and F814W magnitudes (hereafter 777555 and mgu, respectively) 
are expressed in the VEGAMAG systerrQ. The CMD contains 
965 sources going down to 771555 - 29.8 and 111^4 = 29.0. We 
performed artificial star tests (ASTs) to assess the quality of 
the photometric data and evaluate their completeness. The latter 
shows a signi ficant improveme nt with respect to the measure- 
ments of Izotov & Thuan (2004). Not only is the 50% complete- 
ness limit improved by nearly 1 mag, but the completeness is 
also maintained high closer to this limit (the completeness drop 
is more "abrupt"). Furthermore, we used the ASTs to determine 
photometric errors that are more realistic than the default out- 
puts by DOLPHOT/ACS and to evaluate the possible biases in 
the magnitudes retrieved. 

We checked whether the stellar fluxes suf fer from fore- 
groun d Galactic extinction. Whereas the survey of Schleg eTet al.l 
d 1998b yields a Galactic contribution of Ay = 0.11 to 
the extinction toward I Zwl8, the lower-resolution data of 
Burs tein & Heilesl d!982l) indicate an extinction as small as A y — 
0.01. Furth ermore, an analysis of the nebular Ha/H/3 line car- 
ried out by ICannon et al.l d2002l) shows that some regions of the 
main component of I Zwl8 suffer very small extinction, if any. 
Consequently, we assumed that the Galactic extinction toward 
I Zwl8C may be negligible and we did not correct our data for 
it. 



3.2. Selection of the stars 

Before carrying out the analysis of I Zwl8C, we selected the 
stars for use. First, we limited the selection to the stars sit- 
uated above the 90% completeness level. This threshold was 
chosen for two reasons: (i) the uncertainties in the photomet- 
ric analysis (systematic and random errors in the magnitudes, 
evaluation of the completeness map) are small for the region 
of the CMD selected, and (ii) most of the stars situated be- 
low the 90% completeness limit are main-sequence (MS) ob- 
jects that yield little useful information about the age and dis- 
tance of the galaxy, compared with post-MS stars. Then, we re- 
jected the sources situated within the two clusters annotated as 
"C" and "NW" bv llzotov & Thuanl d2004l) . Finally, we removed 
th e five reddest remainin g stars from the selection. According 
to Marigo & Girardj] d2007l) . who revised the evolutionary mod- 
els of asymptotic giant branch (AGB) stars, the presence of very 
red stars in the CMD is possible. However, the characteristics 
of these objects are still very sensitive to the physics used to 
model them. In principle, this problem concerns only extremely 
red stars. As a consequence, we decided to discard the stars with 
colors 777555 - 77i8i4 > 2. The final selection contains 408 stars, of 
which 315 are located in the blue (777555-777814 < 0.6) branch and 
93 in the red one. 

Given the selection procedure, the completeness function p c 
for the selected stars can be written as a function of the "real" 
completeness p* evaluated for the whole CMD: 



See the HST/ACS zeropoints at p c = -j o th g r ^ e °' 9 (12) 

|http : //www . stscF . edu/hst/acs/analysis/ zeropoints] o iuuinl. 
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3.3. Isoch rones 

An important assumption of this work regards the isochrones 
that we decided to use. We adopted the evolutionary tracks of 
iGirardi et all (l2000h available for Z = 0.0004 and computed the 
theoretical absolute magnitudes through the AC S F555W and 
F814W filters^ using the BaSeL 3.1 stellar librarv lWestera et all 
(120021) . 



4. A first model 

4. 1 . Assumptions and procedure 

We proceeded to a first measurement of the distance, extinction, 
and SFH of I Zwl8C. The main assumptions of the calculus were 
the following: 

- a Salpeter IMF (<f>(M) cc AT 2 35 ) in the 0. 15-120 M range; 

- the extinction to be uniform over the stellar population. 
Furthermore, since little work has b een done regardin g the 
extinction in I Zwl8C (see, however llzotov et aLll200ll) . we 
decided to use the extinction Ay as a free parameter; 

- no constraints on the SFH; 

- no concern for the possible effects of binary stars. 

We explored the values of < K pop in the (Am, Ay) space. 
For each examined point of this space, we computed the 
LFs Lj(T;Am,Ay). Using these coefficients, we evaluated 
the corresponding dis tribution N(t) with the algorithm of 
iHernandez et alj ([1999) and derived the entropy < K pop (Am, Ay). 
Furthermore, we calculated the entropies < H& sl (Am) = 
logL dist (Am) and < H ext (A v ) = logL e xt(Ay), where L d ist(Am) and 
£ext(Ay) are the LFs associated to the distance modulus and the 
amount of extinction, respectively. Those two LFs are defined by 

L dist (Am) = I L pop (Am,A v )dA v , (13) 



L eA (A v )= \ Lpo P (Am, Ay) dAm. (14) 



4.2. Results 

The results are presented as follows. In Fig. [3] we present the 
curves of ( H^ l {Am) and ^H^iAy). Figure [4] shows the observed 
CMD along with some theoretical isochrones for the best-fit dis- 
tance modulus and extinction. In Fig. [5] we show the mass rate 
A4(t) derived from N(t) at the best-fit distance and extinction. 

The best-fit distance modulus is Am - 32.12 (corresponding 
to a distance of 26.5 Mpc), and the 99.9% confidence interval for 
Am is 31.71-32.24. Such a distance modulus is higher than any 
previous estimate, by > 1 mag for most of them and a; 0.8 mag 
for the newest (3 1 .4. IXloisi et al.ll2007l) . As for the extinction, 
the best-fit coefficient is Ay = 0.38, and the 99.9% confidence 
interval is 0.30-0.45, a higher value than the different estimates 
of the foreground Galactic extinction (see Sect. 13. 11 1. Thi s value 
falls within the 0.20-0.65 interval derived bv llzotov etaf] d2001l) 
from a spectroscopic study of I Zwl8C and suggests the exis- 
tence of dust clouds responsible for extinction inside this object. 

The estimate of A4(t) shows various interesting features: 



2 Tables of the ACS filter throughputs are available at 
http : //acs .pha . jhu.edu/instrument/photometry 



First model 

_248 r ' ' 1 ' ' ' 1 ' 1 ^ 




-260 t , , ^ i , , , , , J_ ,_ 

31.6 31.8 32.0 32.2 

Am 

First model 
-248 1 I ! 




-258 f .... i .... i .. I . i .... i .... 1 

0.25 0.30 0.35 0.40 0.45 0.50 
Am 

Fig. 3. 'Hdist (upper panel) and "K e xt (lower panel) curves for the 
model of Sect. |4] The dashed lines indicate the best-fit values 
of Am and Ay, while the dotted lines delineate the 99.9% confi- 
dence intervals. 

Firs! model 




Fig. 4. Observed CMD and theoretical isochrones for the model 
of Sect. |4] The observed data are corrected for interstellar ex- 
tinction, using the best-fit value Ay = 0.38. The dotted line 
shows the threshold above which we selected the stars to work 
with. The vertical axes show both apparent and absolute F814W 
magnitude scales, adopting the best-fit distance modulus Am = 
32.12. 
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Fig. 5. SFH derived from the approached explained in Sect. |4] 
The curve was slightly smoothed to limit the apparent statistic 
noise. 



an intense episode of SF is visible at a ges t < 15 Myr. Thi s 
episode has already been detected by llzotov et al.l (12001 1) . 



who analyzed optical s pectra of I Zwl8C, and confirmed by 
llzotov & Thuanl d2004t) in their study of the resolved stellar 
population; 

- a lesser starburst is seen at r ~ 100 Myr; 

- approximately in the 15-70 Myr range the SFR is compara- 
tively moderate but not negligible between those two bursts. 
Interestingly, the corresponding isochrones are those that fall 
between the two branches of the CMD, although they also 
cover populated regions of the CMD; that is, this empty zone 
of the CMD results not from an interruption in the SFH, but 
from the random filling of the CMD. Considering the esti- 
mated SFR, an average of 4 stars would be expected between 
the branches. The likelihood of observing no star there is 
about 6%, which is not small enough to discard the model; 

- the estimate of M(t) suggests a low-rate, continuous SF that 
began at least a few hundreds Myr ago. However, the low 
number of stars benefit from the 300 Myr isochrone, and 
the uncertainty in their ages makes it impossible to confirm 
whether this SF process did occur. 

Although this method seemed to yield accurate measure- 
ments of the distance, extinction, and SFH of I Zwl8C, it was 
important to assess its reliability. Because of the complexity of 
the process, this could only be achieved by means of Monte- 
Carlo simulations. In the next section, we describe how we car- 
ried out this test. 



4.3. Comparison to simulations 

We ran a series of Monte-Carlo simulations of stellar popula- 
tions, using the same assumptions as in Sect. 14.11 and the es- 
timates of the distance, extinction, and SFH described in Sect. 
14.21 For each simulation, we drew 408 stars with random ages 
and masses (using the SFH and the IMF as distribution func- 
tions) and photometric completeness above the 90% threshold. 

We first analyzed 3 simulated populations using a large grid 
of points in the (Am, Ay) space to check the general behavior of 
our algorithm. The three simulations yielded very similar results. 
In Fig.|6j we present one of the synthetic CMDs, Fig.|7]shows the 
corresponding T^ist and 'H ex t curves, and Fig. [8] shows the mass 
rate computed at the best-fit distance and extinction. There is 
no significant discrepancy between the distance, extinction and 




- -s 



m S53 in„ 



Fig. 6. Synthetic CMD and theoretical isochrones for one of the 
simulations Sect. 14. 31 The stars situated below the dotted line be- 
long to the original observations of I Zwl8C, whereas the others 
were simulated. 

SFH estimates obtained for I Zwl8C and for the 3 simulations. 
In principle, this indicates that the probabilistic method we used 
suffers from no severe bias. 

To check for small biases in our approach and verify the con- 
fidence level of our first model, we performed 100 new simula- 
tions, this time exploring a limited area of the (Am, Ay) space 
around the input values. The histograms of 'Hp p,Am - 32.12 
and Ay - 0.38 are shown in Figs. |9l and [TOl They show a possi- 
ble bias in the estimation of Am and Ay, the first quantity being 
slightly overestimated and the second underestimated. However, 
the biases encountered are small (~ 0.03 mag) and do not exceed 
the respective dispersions of the simulations. As for the "Kpop 
distribution, the value obtained for I Zwl8C falls well within the 
range of simulated values, an indication that the observations are 
well-fitted by the model shown in this section. 

Although we have been able to satisfactorily describe the 
CMD of I Zwl8C, it is important to review the assumptions so 
far. It is possible that modifying or relaxing some of them results 
in other satisfactory models of this stellar population. We discuss 
these modifications in the next section. 



5. Towards a better model for I Zw18C 

There are at least four possible sources of error in our modeling 
of I Zwl8C: the isochrones (via the evolutionary tracks, and the 
model atmospheres), the IMF, the treatment of the extinction and 
overlooking binarity. 



5. 1. Effects of the isochrones/stellar luminosity functions 



As explained in Sect. I2.21 our method is related not just with 
the isochrones used, but also with the stellar luminosity function 
they produce: not only is the shape of the isochrone relevant, also 
the IMF (for MS stars) and the lifetime of different evolutionary 
phases (for post-MC stars In this way, the method provides 



3 Note that lifetimes of fast evolutionary phases and the amount of 
stars in a given evolutionary phase are linked by the fuel consump- 
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Fig. 7. 'Kdist (upper panel) and < K e xt (lower panel) curves for one 
of the simulations in Sect. [43] The dashed line indicates the best- 
fit value of Am, while the dotted lines delineate the 99.9% con- 
fidence interval. 

Simulation 



\ 

S 

73 

o 



40 














0.6 












0.5 








30 




0.4 












0.3 












0.2 








20 




0.1 












0.0 












8.5 


9.0 9.5 10.0 j 


10 





















7 


8 




9 


10 






log 


r (a) 





Fig. 8. SFH derived for one of the simulations in Sect l4.3l (bold 
curve). The thin curve represents the SHF estimated for I Zwl8C 
in Sect. 2] and is shown for comparison. 
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Fig. 9. Histogram of 'Hpop for 100 simulations as described in 
Sect. 14.31 The superimposed curve is Gaussian and has the same 
mean and standard deviation as the distribution of <H, 
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Fig. 10. Histograms of Am - 31.12 and A v - 0.38 for 100 simu- 
lations as described in Sect. 14.31 



additional constraint to the ages and distances obtained, as we 
see in Sect. [6] In this situation, just the shape of the isochrones 
from different authors is not enough to assess the uncertainties 
of the isochrones used. 

tion theorem (Renzini & Buzzonilfl984l. 1 1 986L iBuzzonill 1989T) : see also 
iMarigo & Girardil d200lh and lCervino & Lu ridiana ( 2005) for the issues 
about obtain isochrones from evolutionary tracks and its relation with 
the stellar luminosity function. 



However, we can obtain a broad idea of how the results 
would differ by comparing just the mean value of the stellar 
luminosity function produced by different isochrones. By com- 
paring of the results of evolutionary synthesis models (whose 
principal outpu t is actually the mean value of the stellar lumi- 
nosity function, Cervifio & Luridiana 2006). Such a comparison 
has been done bv lBuzzonil (120051) . who in Fig. 1, shows the com- 
parison of the evolution of U, V, K and bolometric luminosity 
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Fig. 11. Function W(t) as described in Sect. 



for different synthesis models (which use different evolutionary 
tracks). In general, there is very good agreement between the 
different synthesis models, hence evolutionary tracks, except in 
the K band (which strongly depends on the synthesis model im- 
plementation of AGB, and post-AGB stars). For the given photo- 
metric bands we used, it is expected that using a different set of 
isochrones produces (at least in a global average) similar results. 

5.2. Effects of the IMF 

Changing the IMF shape modifies the values of the LFs Lj 
through Eq. [2] At a given age, this will be important if the min- 
imum and maximum initial masses "visible" in the photometric 
data, respectively M m m(j) and M lmx {T), are significantly differ- 
ent from each other. Let us examine the case where the IMF is 
a power law of slope a. The ratio between the LFs Lj(j) of a 
■MmaxCr) mass star and a A1 m i n (T) one is nearly proportional to 
(AlmaxM/jVlmtoCT))"- As a result, we can evaluate the impact of 
an uncertainty of Aa in the IMF slope on < Hj(j) with the follow- 
ing function: 



W(t) = Aa log 



M max (T) 
Mnfa(T) 



(15) 



Even though the Salpeter IMF (a = -2.35) is widely em- 
ployed in the study of stellar populations, various authors have 
sugge s ted th e use of other IMFs. For example, according to 
IScalol dl998l) . the IMF slope may vary from a = -2.2 to 
a = -2.8 in the M > 1 M range. Thus, we computed W(t) us- 
ing Aa = 0.5 to check the importance of our choice of a Salpeter 
IMF. The curve is shown in Fig. QT| As seen, at any age (in par- 
ticular at t > 10 Myr), the value of W(t) is low. This reflects 
that our star sample is almost exclusively composed of post-MS 
objects whose properties, at a given age, are very sensitive to 
their initial masses, so the interval of initial masses they cover at 
this age is necessarily narrow. The consequence is that the exact 
choice of the IMF has no significant consequence on the calcu- 
lation of the entropies < Hj(r), and changing the IMF will have 
little impact on the model described in Sect. [4] 

5.3. Effects of binary stars 

Our results may be affected if a significant fraction of the point 
sources used are actually binary stars (either physical or opti- 
cal) where the luminosities of both components are comparable. 



Assume that a given point source is a binary of integrated mag- 
nitude fft8i4, with the magnitudes of the primary and secondary 
component mgi 4j andmgi4,2, respectively (mgu,2 > »Z8i4,i)- The 
magnitude of the secondary companion and the "magnitude ex- 
cess" <5mgi4 = W814 l _ ' w 8i4 caused by its presence are related 
by 



mu,2 = ftgu - 2.5 log{\ - 10 



-QA5m%\4 



)• 



(16) 



The magnitude excess cannot be greater than 2.5 log 2 = 0.753 
because we impose the condition that the secondary companion 
be fainter than the primary one. 

In spite of various attempts, it has not been clearly demon- 
strated that the component masses in physical binaries are cor- 
related (Krou paet al.lll993l) . hence, we assumed that in a given 
binary star, the component masses are independently drawn from 
the IMF. For each point source, we computed the probability Pbin 
of the magnitude excess as Sm^ > <t%u, where <xgi4 is the pho- 
tometric error on m%\n.. We calculated this probability assuming 
the source to be either a physical binary (the two stars have the 
same age) or an optical double star (the ages of the two com- 
ponents are independent). For all the sources of our sample, we 
found Pbin < 0.15%, which means that binarity is very unlikely 
to affect our results seriously. This is actually consistent with 
only looking at the very brightest stars of I Zwl8C: the chance 
for a random star of the galaxy to be nearly as bright as them is 
small. 



5.4. Effects of spacial variations in the intrinsic extinction 

An arguable simplification of our first model is that we assumed 
the interstellar extinction is uniform. In principle, this is not the 
case since dust clouds may be present inside I Zwl8C. The pres- 
ence of such clouds is supported by the estimate of the extinction 
in our first model, and can cause the extinction to vary along the 
projected view of the galaxy. Large-scale fluctuations may ex- 
ist, depending on the morphology of the dust clouds. Moreover, 
even on local scales, the extinction may be distributed, given that 
the stars can be located at different depths inside the clouds. 

To our knowledge, the only work in which a het erogeneous 
extinc tion has been measured in I Zwl8C is that of llzotov et al.l 
(2001). They find the value of Ay to range between 0.20 and 
0.65. This supports the idea that treating the extinction with 
some detail would be useful for more exhaustive modeling of 
IZwl8C. 



6. A model with detailed treatment of the extinction 

We carried out a second model of the stellar content of I Zwl8C, 
based on the same assumptions and procedure as the first model, 
but chaning the treatment of extinction which was dealt with 
in the following fashion. We divided the galaxy into 4 regions, 
shown in Fig. [12] that contain 101, 121, 92, and 59 stars of the 
original selection, respectively. In each region k, we supposed 
the extinction of a given star is a random variable following a 
uniform probability law, with Ay falling in the range AV + AAy\ 
Mathematically, this was achieved by replacing Lj(t; Z,d,A) by 



L'j(r,Z,Am, Af,AAf) 



2AA 



f*+ * L'!{r;Z,Am,A v )6A v 



(17) 



We determined the best-fit values of Ay, AAy\ and Am in 
the following way. For given values of AAL , we computed the 
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Fig. 13. "Kdist curve for the model of Sect. [6] The dashed line 
indicates the best-fit value of Am, while the dotted lines delineate 
the 99.9% confidence interval. 



Second model 



Fig. 12. HST/ACS image of I Zwl8C through the F555W filter, 
with the boundaries of the 4 regions used in Sect.[6]overplotted. 
The two unresolved clusters are also reported. 



Table 1. Best-fit values of AA* and Ay' for the three regions 
described in Sect. [6] 



«) 



Region 




r<« 
A v 


(99.9% CL) a 


1 


0.10 


0.42 


(0.25-0.50) 


2 


0.00 


0.34 


(0.08-0.44) 


3 


0.15 


0.32 


(0.15-0.51) 


4 


0.00 


0.45 


(0.10-0.68) 



The values between parenthesis give the limits of the 99.9% confi- 
dence intervals for A?'. 



fffjAm) curves, derived the <H dist (Am) = % k 'H™ t (Am) distri- 
bution, and picked the best-fit distance modulus Am. At this dis- 
tance modulus, we computed the maximum entropies 1ri {k) as a 
function of Ay. We repeated this operation for different sets of 
AA ( , k ? so as to maximize the values of < H {k) . 



6.1. Results 



The results are presented as follows. In Fig. Qj] we show the 
curve H^P (Am). The best-fit distance modulus is Am = 32.16 
(coiTesponding to a distance of 27.0Mpc) and the 99.9%confi- 
dence interval is 3 1.78-32.29. We found that the best-fit value of 
Am is virtually insensitive to the adopted values of Ay is ac- 
tually close to that found with our first model. In Table Q] we 



summarize the values of A^y 



and AAyK All four extinction aver- 



ages Af fall within the 0.20-0.65 range of llzotov et all (120011) . 
Furthermore, the extinction spreads AAy represent significant 
fractions of the averages for regions 1 and 3. 




1 2 

m 6SS~ m SH 



Fig. 14. Dereddened CMDs for (a) region 1, (b) region 2, (c) 
region 3, and (d) region 4 (see Sect. [6|. The isochrones are the 
same as in Fig. [4] 
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Fig. 15. SFH computed for region 1 (dashed line), regions 2 and 
3 (dotted line), and all three regions (full line). 
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Fig. 16. Position of the stars falling in the following age ranges: (a) 015 Myr, (b) 1560 Myr, (c) 60120 Myr. For reference, a surface 
brightness isocontour of the galaxy was also plotted in the panels. 



6.2. The SFH of I Zw18C and its spatial distribution 

In Fig. [TU we show the CMDs for the individual regions at the 
adopted distance and average extinctions. Figure [15] shows the 
mass rate curves for the four regions together. Finally, Fig. Q~6] 
shows the position of the stars in different age ranges. The age 
Tj of star j was defined as the age maximizing the LF L"(t). 

The global SFH shown in Fig.[l5]is nearly identical the one 
discussed in Sect. [4] In particular, the two starbursts at r < 15 
Myr and r ~ 100 Myr and the continuous SF between them 
are confirmed with our new model. However, it is evident that 
the SFH is not homogeneous across the galaxy. The young stars 
(< 15 Myr) tend to concentrate in the southwest tip of the ob- 
ject, while there is a slight excess of 60-120 Myr old stars in the 
northeast lobe. On the other hand, intermediate age stars (15- 
60 Myr) are distributed rather evenly through the component of 
the galaxy. This observation suggests that I Zwl8C forms stars 
both through a continuous, global process and through local star- 
bursts. The existence remains unconfirmed of the longlasting SF 
that started a few hundreds to thousands of Myr. Its possible con- 
tribution to the total mass of the stellar population is very small, 
about 10~ 3 , and I Zwl8C can be considered mostly as a young 
object. 

6.3. The age of I Zw18C 

So far, we have presented estimates of the SFH of I Zwl8C 
where the stars could be as old as 14 Gyr. These estimates do 
not provide a precise constraint on the age of the stellar sample 
considered, i.e. the age of the oldest of its stars. To establish the 
age, we evaluated the distribution 7V(t) for a given tested age 
Tpop, with the coefficients Gj(j) defined as 

Gj(r) = j j P{Am,Af)P'^T;hm,Af)6Am6Af, (18) 

forcing N(t) — at r > T pop . We then retained the corresponding 
entropy < H a ge(T pop ) and derived the LF L age (T pop ). The definition 
of Gj{t) has the advantage of properly accounting for the proba- 
bilistic distributions of Am and Ay\ With a few simulations, we 
checked that this procedure yields the correct LF. 

In Fig. [17] we show the curve L- dge (T pop ) obtained. If all the 
stars are considered as valid elements of the sample, then the 
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Fig. 17. Likelihood function of the age of I Zwl8C, as computed 
in Sect. 1631 

most likely age of I Zwl8C is 125 Myr, with a 99.9% confidence 
level interval of 1 10 Myr-7.6 Gyr. To ensure that (r pop ) < 1 Gyr, 
it would be necessary to remove the oldest 20 stars from the 
sample. That is, I Zwl8C is very likely a young object, but we 
cannot rule out the possibility that this object is several Gyr old. 

6.4. A high distance estimate 

Our estimate of the distance to I Zwl8 is by far greater than 
those assumed or evaluated in all previous works published on 
this object, even considering the error bars of the different esti- 
mates. No modification to our model, including the removal of 
a large number of stars from the original sample, could alleviate 
this discrepancy. Whatever stars were removed, the distance es- 
timate was found to be very similar to that described in Sect. 16. II 
A good description of the distribution of the colors and magni- 
tudes of the stars can be obtained by assuming the distance to 
I Zwl8C to be 27 Mpc, as shown in Fig. [18] In particular, the 
fit of the magnitude distribution is excellent. Let us now adopt 
the distance modulus o f 31.3, inferred from the RBG population 
by lAloisi et al.l d2007l) . As shown in Fig. [19] where the extinc- 
tion and the SFH were re-fitted, a significant excess of luminous 
AGB stars (772555 - '«8i4 ~ L7, 111^4 « 24) is predicted, while 
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Fig. 18. Color and magnitude histograms of the stars considered 
in Sect. [6] The thick curves represent the average distributions 
expected at the best- fit distance and extinction coefficients. 
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Fig. 19. Same as Fig.Q~8]but at a distance modulus of 31.3. 



the expected number of MS stars at the faint end of the CMD 
(m555-msi4 ~ 0, ;«8i4 ~ 28) is lower than the observed one. This 
scheme worsens if we assume an even shorter distance. Thus, it 
appears that the large distance modulus we computed is essen- 
tial for explaining the distribution of colors and magnitudes of 
the stars observed. 



7. Conclusions 

In an attempt to evaluate the distance to the dwarf galaxy 
I Zwl8C and its SFH, we carried out a detailed analysis of its 
resolved stellar content. We developed a probabilistic approach 
based on the current knowledge about SF and its random aspects, 
in order to exploit as much of the information held by the pho- 
tometric data as possible. We used high-quality HST/ACS ob- 
servations of the galaxy, on which we performed state-of-the art 
photometry to obtain high signal-to-noise and high completeness 
data. 

We carried out a first model with a few simple assumptions 
about the stellar population studied. We found the distance to 
I Zwl8C to be as high as 27 Mpc and detected several impor- 
tant features in the SFH: two starbursts of which an ongoing one 
(r < 15 Ma) dominates the SFH and the other occurred a; 100 
Myr ago, a rather constant SF episode that has lasted for 100 
Myr or so, and a possible low-rate SF process for larger ages. 
We performed Monte-Carlo simulations to assess those results. 
Although the simulations proved our approach to be efficient, 
some of our assumptions could limit the reliability of the model. 
We reviewed the assumptions and found that only the treatment 



of the extinction could yield a more robust depiction of the stel- 
lar population studied. 

We performed a second model, accounting for the spatial 
variations and random fluctuations of the interstellar extinction 
inside I Zwl8C.We could confirm the high estimate of the dis- 
tance modulus of I Zwl8C (Am = 32.16+°^, i.e. d = 27.0+^ 
Mpc), significantly greater than any value previously published. 
The description of the SFH was also maintained. The spatial dis- 
tribution of the stars in various age intervals suggests that the 
stellar content of I Zwl8C has formed both through a global, 
continuous process and through a series of local starbursts. 

That the distance we computed for I Zwl8C is greater than 
those published in previous works is important for study the re- 
solved stellar population of both this object and the main body of 
I Zwl8. For example, the spectra assigned to the ionizing stars 
will be more intense if one adopts our distance instead of a more 
conservative one. The age estimate of the red stars is also sen- 
sitive to the adopted distance, because it is more at shorter dis- 
tances. Anothe r consequence of our d istance estimate is that it 
explains why llzotov & Thuanl (2004) did not detect any RGB 
stars in the galaxy. 

An equally important result of our work regards the age of 
the galaxy. Our probabilistic approach allowed us to state that 
I Zwl8C is most likely young, with a maximum likelihood age 
of 125 Myr, although we could not discard ages of several Gyr. A 
more accurate age could be inferred from photometric data that 
is either deeper or that covers a wider range of infrared wave- 
lengths. In any case, our estimate of the SFH shows that nearly 
all the stellar content of I Zwl8C has formed in the last few hun- 
dred Myr. 

Let us finish the paper with a final note about how its re- 
sults have affected the authors of the paper (and maybe also the 
reader). Each of us had own expectations about this work. Some 
of us expected a shorter distance for I Zw 1 8 C and (by ex t ension 
to the main body) more consistent with pre- llzotov et al.l ([200 1 ) 
estimations: a lower distance implies one lower absolute lumi- 
nosities and a higher impact from IMF sampling effects in the 
I Zwl8 system studies. Others expec ted that this study wo uld 
produce similar results to the ones bv llzotov & Thuanl (2004). 

In all cases, our preconceptions (based in different readings 
and interpretations of previous works and our different personal 
history in I Zwl8 studies) have not been satisfied. 

Among us, a posted criticism has arisen (maybe without be- 
ing completely aware of our expectations and its impact on the 
interpretation of the results). As an example of potential problem 
that the reader might also have thought about, part of our analy- 
sis is based only in a few luminous red and MS stars. However, 
our solution is based on an overall analysis of both populations 
(i.e. the luminosity function) and not on a few stars in a par- 
ticular box of the CMD in different areas (Figs. [T8l and [191 and 
discussion in Sect. 16.4b . We are not able to find any mistake in 
the methodology or in the results, and any possible criticism also 
applies (more strongly) to results from standard CMD method- 
ology. 

Independently of the particular expectations of each of the 
authors, we hope that the controversy about the I Zwl8 system 
distance and I Zwl8C star formation history will be solved with 
deeper observations of I Zwl8C using future large observational 
facilities. But for the moment, against our expectations, it seems 
that current data point toward a distance of 27 Mpc to the I Zw 1 8 
system and maximum likelihood age of 125 Myr. 

Acknowledgements. This work was supported by the Spanish Programa 
Nacional de Astronomia y Astroflsica through FEDER funding of the project 



12 



L. Jamet et al.: The distance and SFH of I Zw 18 C 



AYA2004-02703 and AYA2007-64712. LJ was supported by a UNAM post- 'I Zwl8C on page[l] 
doctoral grant. 



References 

Aloisi, A., Tosi, M., & Greggio, L., 1999, AJ, 118, 302 
Aloisi, A., et al. 2007, ApJ, 667, L151 

Brown, T. M., Heap, S. R., Hubeny, I., Lanz, T., & Lindler, D. 2002, ApJ, 581, 
L129 

Burstein, A., & Heiles, C, 1982, AJ, 87, 1 165 
Buzzoni, A. 1989, ApJS, 71, 817 
Buzzoni, A. 2005, MNRAS, 361, 725 

Cannon, J. M., Skillmann, E. D., Gamett, D. R., & Dufour, R. J., 2002, ApJ, 565, 
931 

Cervino, M., Luridiana, V., & Castander, F. J. 2000, A&A, 360, L5 
Cervino, M., Valls-Gabaud, D., Luridiana, V., & Mas-Hesse, J. M. 2002, A&A, 
381, 51 

Cervino, M., & Luridiana, V. 2005, in Resolved Stellar Populations, 
D. Valls-Gabaud & M. Chavez (eds.), ASP Conf. Ser. (in press), 
|arXiv:astro- ph/051041 1 

Cervino, M., & Luridiana, V. 2006, A&A, 451, 475 

Dolphin, A., 2000, PASP, 112, 1383 

Dufour, R. J., Esteban, C, & Castaneda, H. O., 1996a, ApJ, 471, L87 
Dufour, R. J., Garnett, D. R., Skillman, E. D., & Shields. G. A., 1996b, ASPC, 
98, 358 

Gallart, C, Freedman, W. L., Aparicio, A., Bertelli, G, & Chiosi, C. 1999, AJ, 
118, 2245 

Girardi, L., Bressan, A., Bertelli, G, & Chiosi, C, 2000, A&AS, 141, 371 

Hensler, G, & Rieschick, A., 2002, ASPC, 285, 341 

Hernandez, X., Valls-Gabaud, D., & Gilmore, G. 1999, MNRAS, 304, 705 

Hunter, D. A., Elmegreen, B. G, & Baker, A. L., 1998, ApJ, 493, 595 

Izotov, Y. I., Chaffee, F. H., Foltz, C. B., Green, R. R, Guseva, N. G.,& Thuan, 

T. X. 1999, ApJ, 527, 757 
Izotov, Y. I., Chaffee, F. H., Foltz, C. B., Thuan, T. X., Green, R. F, Papaderos, 

P., Fricke, K. J., & Guseva, N. G. 2001, ApJ, 560, 222 
Izotov, Y. I., Foltz, C. B., Green, R. F, Guseva, N. G, & Thuan, T. X. 1997, ApJ, 

487, L37 

Izotov, Y. I., Stasihska, G, Guseva, N. G, & Thuan, T. X. 2004, A&A, 415, 87 

Izotov, Y. I., Thuan, T. X., & Lipovetsky, V., 1997, ApJS, 108, 1 

Izotov, Y. I., Thuan, T. X., & Stasihska, G. 2007, ApJ, 662, 15 

Izotov, Y. I., & Thuan, T. X. 1998, ApJ, 497, 227 

Izotov, Y. I., & Thuan, T. X. 1999, ApJ, 511, 639 

Izotov, Y. I., & Thuan, T. X. 2004, ApJ, 616, 768 

Kroupa, P., Tout, C. A., & Gilmore, G. 1993, MNRAS, 262, 545 

Kunth, D., & Ostlin, G, 2000, A&ARev., 10, 1 

Legrand, F, Kunth, D., Roy, J.-R., Mas-Hesse, J. M., & Walsh, J. R. 1997, A&A, 
326, L17 

Legrand, F, Kunth, D., Roy, J. R., Mas-Hesse, J. M., & Walsh, J. R., 2000, A&A, 
355, 891 

Marigo, P., & Girardi, L. 2001, A&A, 377, 132 
Marigo, P., & Girardi, L. 2007, A&A, 469, 239 
Meynet, G, & Maeder, A. 2005, A&A, 429, 581 
Ostlin, G, 2000, ApJ, 535, L99 

Pagel, B. E. J., Simonson, E. A., Terlevich, R. J., & Edmunds, M. G. 1992, 
MNRAS, 255, 325 

Papaderos, P., Izotov, Y. I., Thuan, T. X., Noeske, K. G, Fricke, K. J., Guseva, 

N.G. & Green, R.F. 2002, A&A, 3939, 461 
Peimbert, M., & Torres-Peimbert, S. 1974, ApJ, 193, 327 
Peimbert, M., & Torres-Peimbert, S. 1976, ApJ, 203, 581 
Peimbert, M., Luridiana, V., & Peimbert, A. 2007, ApJ, 666, 636 
Renzini, A., & Buzzoni, A. 1984, Spectral Evolution of Galaxies, 86 
Renzini, A., & Buzzoni, A. 1986, Spectral Evolution of Galaxies, 122, 195 
Robertson, B., Bullock, J. S., Font, A. S., Johnston, K. V., & Hernquist, L., 2005, 

ApJ, 632, 872 

Scalo, J., 1998, The stellar initial mass function, ASP Conf. Ser., 142, 202 
Schlegel, D.J., Finkbeiner, D.P. & Davis, M., 1998, ApJ, 500, 525 
Skillman, E. D., & Kennicutt, R. C, Jr. 1993, ApJ, 411, 655 
Tosi, M., Aloisi, A., Mack, J., & Maio, M. 2007, IAU Symposium, 235, 65 
van Zee, L., Westpfahl, D., Haynes, M. P., & Salzer, J. J., 1998, AJ, 1 15, 1000 
Vilchez, J. M., & Iglesias-Paramo, J. 1998, ApJ, 508, 248 
Westera, P., Lejeune, T, Buser, R., Cuisinier, F, & Bruzual, G. 2002, A&A381, 
524 



List of Objects 

'I Zwl8' on pageQ] 



